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The one-loop induced top quark decay t —f u\U2U2 {ui — u,c) is calculated in the context of 
the standard model. The dominant contribution to this top quark decay arises from the Feynman 
diagrams induced by the off-shell tuig* vertex, whereas the box diagrams are negligibly small. In 
contrast with the on-shell tuig vertex, which only gives rise to a pure dipolar effect, the off-shell 
tuig* coupling also involves a monopolar term, which gives a larger contribution than the dipolar 
one. It is found that the branching ratio for the three-body decay t U1U2U2 is about of the same 
order of magnitude of the two-body decay t —> uig, which stems from the fact that the three-body 
, decay is dominated by the monopolar term. 
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; I. INTRODUCTION 

The top quark detection at the Fermilab Tevatron |Q greatly boosted the interest in top quark physics. The 
^ • large mass of this quark suggests that it could be very sensitive to new physics effects, which may manifest themselves 
■ through anomalous rates for its production and decay modes. Although some properties of the top quark have already 
t been examined at the Tevatron 0, a further scrutiny is expected at the CERN large hadron collider (LHC). This 
^ machine will operate as a veritable top quark factory, producing about eight millions of tt events per year in its first 
00 . stage, and hopefully up to about eighteen millions in subsequent years Yet in the first stage of the LHC, many 
rare processes involving the top quark are expected to be accessible. It is thus worth investigating all of the top quark 
[ decays within the standard model (SM) in order to find out any scenario that may be highly sensitive to new physics 
^ ■ effects. 

I In the SM, the main decay channel of the top quark is t bW. Although the nondiagonal t — > dW and t sW 
. modes are more suppressed, they still have sizable branching ratios. For instance, Br{t sW) is of the order of 
10^'^. As far as rare decays are concerned, the three-body tree-level induced modes t diWZ and t — > uiWW, 
with di — b,s,d and ui — u, c, are strongly dependent on the precise value of the top quark mass. It has been shown 
that the t uiWW decays are severely GIM-suppressed 0, but t bWZ can have a branching ratio of the order 
Qh! of 10~^ for a top quark mass larger than 187 GeV 0. This decay mode has been suggested as a probe for the top 
^ ' quark mass because it is almost in the threshold region At the one- loop level, there arise the flavor changing 
neutral current (FCNG) decays t uiV {V = QjlZ) and t — > uiH, which are considerably GIM-suppressed, with 
branching ratios ranging from 10^^" to 10^^"^ |3iB>l3- Motivated by the fact that any process that is forbidden or 
strongly suppressed within the SM constitutes a natural laboratory to search for any new physics effects, FCNG top 
, quark decays have been the subject of considerable interest in the literature 1 0., Jd, .12 , ,14j . It turns out that they 
d may have large branchin g ra tios, much larger than the SM ones, within some extended theories such as the two-Higgs 
doublet model (THDM) |lOl| , supersymmetry (SUSY) models with nonuniversal soft breaking , SUSY models with 
broken i?-parity [l^ . and even more exotic scenarios ^13]. Similar results for the decays t — > uiV and t — > uiH were 
obtained within the context of effective theories p^ . 

In this work, we present a calculation of the t — > U1U2U2 decay {u2 stands for the u or c quark), which arises at the 
one-loop level in the SM. Although the study of rare top quark transitions has attracted considerable attention, to our 
knowledge the rare decay t —^ U1U2U2 has never been analyzed before. The rest of the paper is organized as follows. 
Section II is devoted to the analytical calculation of the decay t — > M1U2U2. The numerical results and discussion are 
presented in Sec. HI along with the final remarks. 



II. THE DECAY t U1U2U2 



Decays of the type t — > U1U2U2 proceed through the reducible diagrams shown in Fig. ^i), which are mediated 
by the Z, H, ^ or g bosons. While those Feynman diagrams mediated by the Z and H bosons are enhanced due 
to the fact that the intermediary boson is on resonance (provided that niH < Tit), those diagrams mediated by the 
photon (gluon) are enhanced by the effect of the photon (gluon) pole. There are also contributions from box diagrams 
carrying W gauge bosons and down quarks [see Fig. ^ii)]. Since each type of diagram renders a finite amplitude 
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FIG. 1; Feynman diagrams contributing to the t u\U2U2 decay. The bubble stands for all the contributions of the type 
shown in Fig. |21 
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FIG. 2: Feynman diagrams contributing to the tuig vertex in the unitary gauge. In the -Rj gauge there are an extra set of 
diagrams in which the W boson is replaced by its associate pseudogoldstone boson. 

by its own, the different contributions can be considered as independent. It can be shown that the t — s- U1U2U2 
decay is essentially determined by those graphs involving a virtual gluon, i.e., those reducible diagrams involving the 
one-loop vertex tuig* and the tree-level vertex g*U2U2- The role played by this contribution is evident from the fact 
that it constitutes an electroweak-QCD mixed effect. This is to be contrasted with those reducible graphs mediated 
by the 7, Z and H bosons, as well as the box diagrams, which are entirely determined by electroweak couplings. 
As a consequence, the pure electroweak contributions become suppressed by a factor of a/ as as compared with the 
electroweak-QCD mixed ones. 

Once the most relevant properties of the t — > U1U2U2 decay were described, we would like to emphasize a noteworthy 
feature of this process. It is closely related to the fact that the t U1U2U2 decay is mediated by a virtual massless 
vector boson, i.e., the gluon or the photon. Without losing generality, it is enough to discuss the gluon contribution 
as it is the dominant one. Naively, one would expect that the rate for the two-body decay t — > uig is larger than 
that for the three-body decay t — > U1U2U2, which in fact is not true. This stems from the fact that while the on-shell 
tuig vertex is characterized by a dipole structure (the tui pair couples to the gluon through the gauge tensor GJ^^^), 
the corresponding off-shell tuig* vertex involves also a monopole structure (the tui pair interacts directly with the 
gauge field). Therefore, while the t uig transition is entirely determined by the dipole structure, both the 
dipole and the monopole structures contribute to the t — > U1U2U2 process. It turns out that the contribution from the 
monopolar term can be considerably larger than that arising from the dipolar one. We have found that this is indeed 
the case for the rare decay t — > uiU2U2. It means that while the t uig decay is determined by the dipolar term, 
the t U1U2U2 mode is governed by the monopolar one. Moreover, the three-body decay is unsuppressed because 
it includes the QCD vertex g*U2U2, which is much less suppressed than the electroweak vertices "f*U2U2 and Z*U2U2. 
The above properties nicely conspire to enhance the t —^ U1U2U2 decay rate by about one order of magnitude with 
respect to that of the t —>■ uig transition. 

We now turn to analyze the general structure of the tuig vertex, which is generated at the one-loop level via the 
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Feynman diagrams shown in Fig. [5] The most general form of this vertex involves up to ten form factors, which 
are associated with the Lorentz structures j^Pl,r, {p - Pi)^iPl,r, q^iPL,R, tT^!/(p - PiYPl,r, and cF^uq" Pl.r, with 
Pl.r = (1 T 15) I'^- There are however a few independent form factors. After imposing the on-shell conditions on 
the fcrmionic fields, the Gordon identities allows one to eliminate the four form factors associated with {p — pi)fiPL.R 
and (Jfii^ip — PiY Pl,r- In addition, the q^e^^ = condition, valid for a real gluon, can be safely used as the off- 
shell gluon couples to a pair of approximately massless quarks. Thus, the only non negligible contributions to the the 
t — > U1U2U2 decay are those of the monopolar {■^fj,PL,R) and dipolar {af^uq'^ Pl,r) terms. Furthermore, no contributions 
proportional to Pr can arise in the = limit. Consequently, the vertex function associated with the tuig coupling 
can be written as 



2 ^ 



Fi{q^Y,,PL + —F2{q')PL<J^,q'' 
nit 



(1) 



where gs is the coupling constant and A°/2 stand for the generators associated with the color group. It is worth 
mentioning that the monopolar contribution Fi vanishes in the on-shell limit as a consequence of the Ward identity 
g^r^ = 0. This means that a FCNC vertex involving an on-shcU gluon can only arise through a dipolar term, such 
as occurs in electrodynamics. This is not true for an off-shell gluon: in such a case the monopolar term yields the 
dominant contribution. This behavior will be explicitly shown below. We will first calculate the form factors in an 
gauge and, in order to assure that our result is gauge-independent, we will also calculate such form factors via the 
unitary gauge. In the gauge the calculation leads to the following amplitude 



where 



2 2^yu,^KtJ (2^)15 A,' 



Tl = 1pPl{^ +1/2 + m,Y^{^ + tti^Y^PlP''^, 
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In the previous expressions, the m„j = approximation was used. In addition, denotes the mass of the internal 
down quark and ^ is the gauge parameter. For simpUcity the calculation was performed in the t'Hooft-Feynman 
gauge = 1). As a crosscheck, we also have performed the calculation via the unitary gauge(^ = oo), in which 
there are only contributions from the first three terms in Eq. (|2Jl, with P"^ replaced by P"^ — g"^ — k°'k^ /rn^. 
The results obtained by these two calculation schemes do coincide, which guarantees that the form factors associated 
with the monopolar and dipolar structures of the tuig* vertex are gauge- independent. Introducing the definition 
Fi = {a^ /^Ti)Ai we can write 

= 2x a-x) ^ V° + ^ ^-^'^"^"^ + ""'^ ' ^^^^ 

^ i—d,s,b \ a—1 / 

M = ^ ^ r K.^^tU /o' + E/'^o(«) + 2™?52Co(a;,a;,) ) , (17) 

ixw{^ X) ^^^^^^ \^ J 

where x = q^/rri^, xw = Tny^/ml^ and xi ~ mf/rri^. The functions depend on Xi and read 

/O = (18) 



and 



/i = j^(^{2{l-x) + ix^A)ix,+xw))x^ 

- 2(2(1 - x) + {x ~ 4:)xw)xw^ , 
fl = Mx + 2)(^1 - X - 2(a;j + xvf))^;^ 



2 



JO ~ 



(19) 



(2(0; + 2)xw - a;(l - a;) - 2(1 - x^)) xw j (20) 



fl = ■p-^-^^(-x2 + 5x-4 + (x^ -3x + 8)(xj+xw))xj 

- 2 (^x^ + 3x - 4 + (x^ - 3x + S)xw)x^^ , (21) 

gi = ^2(x(l -x)'^ + (x + 2)x^ + 2(x^ - l)xw)xw 

+ ((l-x)2(x, -l) + (x + 2)(x2-3a;2,.)-2((x + l)2-4)xH.)x,J, (22) 



(23) 



fl = Y^^(2(l-x)-(x + 2)(x, +xw))x. 

!(2(l-x)-(x + 2)xH/)xiv), (24) 
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+ 2(6xvv-5(1-x))xh/^, (25) 

fi = ^ ^(^x^ + X - 2 - (x^ - 5x - 2){x^ + xwi^Xi 
— 2^3x^ — X — 2 — {x^ — 5x — 2)xw^xw^ , 



(26) 



92 = _ \ 2x[{\ ~ .t)(1 - X - Axw) + Sxl^^xw 

+ (^3.ra,^ + (x, -1)(1-x)^ + (2(1 + 2x)-3x(2.t + 3xh/))xh')x,^, (27) 

In writing the above expressions, the unitarity condition X^i^uii^l ~ was taken into account, i.e., any term 
independent of the internal quark mass was dropped out. Also, it is straightforward to show that J2l=i /a = fo^' 
& = 1, 2, which means that, as expected, the Ai amplitudes are free of ultraviolet divergences. 

We now are ready to calculate the contribution of the tuig* vertex to the t — > U1U2U2 decay. Below, p2 and p2 will 
stand for the 4-momenta associated with the U2 and U2 quarks. It is useful to introduce the following dimensionless 
variable y — {pi +^2)^/™?- The ui quark mass will be retained in the phase space integral since a factor of 
associated with the gluon pole, enters into the t U1U2U2 squared amplitude. Using the expressions for the tuig* 
and g*U2U2 vertices, it is straightforward to construct the invariant amplitude associated with the diagram (i) of Fig. 
^ After making this, we can write the invariant mass distribution dT/dx as follows 

dT{t U1U2U2) mt 



dy V \M\\ (28) 



with the squared amplitude being 



dx 2567r-^ / ^-^ 

spins 



2 2 



^ \M\' - -^(Fi(x,y)|Aip + f^i2(x,y)2i?e(yliA*) 



spins 



|2\ 



+ F2ix,y)\A2\'), (29) 
As far as the Fi(x, y) functions are concerned, they are given by 

Fi(x,y) = -4(x2+x(2y-l) + 22/(2/-l)), (30) 

Fi2(x,2;) = -4x(l-x), (31) 

F2(x, y) = 4x(2y(y - 1) + x{2y - 1) + l) . (32) 

whereas the Passarino-Veltman scalar functions Bq and Cq are, in the usual notation: 

Bo(l) -So(0,m2x„m^), (33) 

i?o(2) = Bo{mfx,m^Xi,mfxi), (34) 

Bo(3) = Bo(TOj,mjXi,TO^), (35) 

Co(x, Xj) = Co(mt , 0, m^x, m^Xi, m^, mfxi). (36) 

The integration limits are as follows 



yrmn = ^0 + ^ (1 " ^) ( 1 " \/ 1 " ^ ) , (37) 
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^yl<x<{l- xo)^ (39) 

where xq — mu-^^/rrit and j/o — ^u-ilfnt- 

From Eq. it is evident that Ai, the monopole term, vanishes in the on-shell hmit (a; — > 0), in agreement 

with the fact that the t — > u\g decay is only determined by a dipole term. We wiU show below that the monopole 
contribution is slightly larger than the dipole one. 

Since the Ai amplitudes do not depend on y, this variable can be integrated over readily. In the ^ limit, we 
are left with 

— - "'"'f * 3 + hi{x)Re{A^Al) + hi.x)\A^\^) , (40) 



dx 1728 s,^ 



where 



fi{x) = ^{l + 2x){l-x)\ (41) 
fi2 = (42) 

X 

f2{x) = -{2 + x)(l-xf. (43) 

X 

It is interesting to note that we will not take into account the limiting case yo ^ {^U2 = 0) when integrating 
over X because the dV /dx distribution would become undefined in a; = due to the gluon pole. This corresponds to 
the case when the U2 quark emerges parallel to U2 and we cannot take the limit of massless U2 quark as it would lead 
to a collinear singularity. Thus, although we have neglected the outgoing quark masses in the transition amplitude, 
they must be retained in the integration limits of the x variable. 



III. NUMERICAL RESULTS AND FINAL REMARKS 



For the numerical analysis we will use the values of the running coupling constant as and quark masses at the rrit 
scale, namely, as{mt) = 0.10683, mt{mt) = 174.3 GeV, mi^imt) — 2.85 GeV, mc{mt) = 0.63 GeV, ms{mt) — 0.09 
GeV, md{mt) — 0.0049 GeV, and TO„(TOt) = 0.00223 GeV 15]. It is worth noting that the numerical results do not 
change considerably for small variations of the outgoing quark masses. 

We first would like to compare the size of the off-shell tuig* vertex with that of the on-shell one. Numerical 
evaluation shows that the tuig dipole contribution is two orders of magnitude smaller than the tuig* monopole 
contribution and one order smaller than the tuig* dipole contribution. Thus, while the t uig decay only receives 
the contribution of the dipolar term, the t — > U1U2U2 transition receives an extra contribution of the monopolar term, 
which is slightly larger than the dipolar contribution. 

The fact that the contribution of the monopole form factor to the tuig* vertex is larger than that of the dipole one 
is exhibited in the invariant mass distribution dT(t — > uiU2U2)/dx, which is shown in Fig. |31 where we have plotted 
separately the monopolar and dipolar contributions. Therefore, it is evident that the t — > U1U2U2 decay is slightly 
dominated by the monopolar term. 

We now turn to the mmrerical evaluation of Br{t uiU2U2). Using the values given above and T{t — > bW) — 1.55 
GeV, we obtain 

Br{t U1U2U2) = 3.38 X 10~^^ (44) 

On the other hand, according to the literature Br{t — > uig) — 5.73 x 10~^^ 0,E|- This result shows that Br{t — > 
U1U2U2) is about of the same order of magnitude than Br(t —> uig). If one sums over all the possible U2U2 pairs, the 
resulting Br{t — > U1U2U2) is of the order of 10"^"'^ and thus larger than BR{t — > uig). Although these decay rates 
seem exceedingly small to be detected ever, they may be largely enhanced in some SM extensions. In such a case the 
effect discussed above may have some interesting implications. 

In conclusion, we have shown the interesting fact that three-body decay t U1U2U2 has a branching ratio about 
the same order of magnitude than the one of the two-body decay t — > uig. Although rare decays of this type are very 
suppressed in the SM, they may have much larger branching ratios in other SM extensions, thereby constituting an 
interesting place to search for any new physics effects. 
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FIG. 3: Monopolar and dipolar contributions to the t U1U2U2 invariant mass distribution dT/dx. 

Note added 

After this work was submitted, a preprint was posted to the preprint archive by Eilam, Frank and Turan[T^. who 
evaluate the t — > egg and t U1U2U2 decays. Although these authors do not presente explicit analytical expressions 
for the t — > U1U2U2 decay |16|. our numerical result for the branching ratio agrees with theirs. We have also learnt 
that Deshpande, Margolis and Trottier presented a similar analysis in These authors reached a similar conclusion 
on the t — > cqq decay in both the standard and the two-Higgs doublet models. 
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